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Abstract 

Background: Abdominal visceral fat affects the metabolic processes, and is an important risk factor for morbidity 
and mortality. The purpose of the study was to develop a quick and accurate estimate in the visceral fat area (VFA) 
of the L4-L5 vertebrae using anthropometric predictor variables that can be measured conveniently. 

Methods: A total of 227 individuals participated in this study and were further divided into a IVlodeling group 
(MG) and a Validation group (VG). Anthropometrics measurements (height, weight, waist circumference, hip 
circumference, age, and subcutaneous fat thickness) and VFAcj were measured using computer assisted 
tomography for all participants. Multivariate linear regression analysis was applied to the MG to construct a 
VFA estimator using anthropometric predictor variables and to evaluate its performance using the VG. 

Results: The estimate equation obtained from the MG were VFAant = -1 44.66 + 1 .84X, + 1 .35X2 + 0.52X3 
(r = 0.92, SEE =14.58 cm^, P< 0.001, n = 152). The Xi, X2, and X3 variables in the equation were denoted as waist 
circumference (WC), age, and abdomen subcutaneous fat thickness (AS). In addition, the correlation between 
VFAant and VFAct showed a high correlation (r = 0.92). 

Conclusion: A rapid and accurate VFA estimation can be achieved by using only age, WC, and AS. The approach in 
the present study provides an easy and reliable estimate that can be applied widely in health and epidemiology 
studies. 

Keywords: Computed tomography. Stepwise linear regression analysis. Subcutaneous fat thickness. 
Waist circumference. Waist-to-hip ratio 



Introduction 

Abdominal fat can be divided into visceral fat and sub- 
cutaneous fat. Previous literatures suggest that the abdom- 
inal visceral fat affects the metabolic processes and is an 
important risk factor for morbidity and mortality [1-3]. 
The abdominal fat can be precisely and reliably measured 
using magnetic resonance imaging (MRI) or computed 
tomography (CT) [4,5]. However, the procedure is not 
only expensive but also limited to hospitals or research 
centers for clinical and epidemiological studies. To over- 
come these problems and reduce costs, current research 
seeks to find more accurate and precise indicators of 
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anthropometric abdominal fat [6-8]. Commonly used an- 
thropometric measures to estimate abdominal visceral fat 
include body mass index (BMI) [9], waist-to-hip ratio 
(WHR) [10], waist-to-thigh ratio [10], waist circumference 
[8,11,12], sagittal abdominal diameter [8,13] and skinfold 
thickness [14]. 

Although those anthropometric measures are easy to 
obtain and are usually highly reproducible, studies that are 
associated with the research in abdominal visceral fat area 
(VFA) uses one or more anthropometric measures are in- 
conclusive. Therefore, by using different anthropometric 
variables to estimate VFA provides different interpreta- 
tions and physiological significance [13,15]. On the other 
hand, when too many anthropometric predictor variables 
are introduced on the right-hand side of a regression 
model, coUinearity becomes an issue. An overly complex 
model, of which there are more predictor variables to 
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collect, will not necessarily improve the accuracy of esti- 
mation [14]. Therefore, it is important to find an optimal 
set of anthropometric measures that would lead to a high 
accuracy rate for VFA estimation yet would be easy to ob- 
tain and to be reliably interpreted. 

In the present study, the CT data on the amount of fat 
around the human L4-L5 vertebrae were used to esti- 
mate VFA (visceral fat area) [1]. Stepwise multivariate 
linear regression was applied to the modeling group to 
build a VFA estimate using BMI, height, age, WHR, WC, 
and subcutaneous fat thickness at several sites as inde- 
pendent predictor variables. The VFA estimation model 
was then tested using data from the validation group. 

Subjects and methods 

Subjects 

A total of 227 Taiwanese adult subjects (Asian) were re- 
cruited voluntarily and agreed to participate by signing 
informed consent forms. A subject pool was generated 
by posting advertisements about the nature of the study 
throughout Taiwan. Participants filled out survey infor- 
mation that included demographic information, their 
physical characteristics and their health status. None of 
the participants had an endocrine, nutritional or growth 
disorder and did not have any major chronic diseases. 
Individuals that had been diagnosed with hypertension, 
diabetes, cancer, renal dysfunctional liver disease or 
long-term pulmonary asthma or who were pregnant 
were excluded from the study. This project was imple- 
mented from the Department of Radiology at Taichung 
Jen-Ai hospital and approved for ethical approval by the 
Institutional Review Board at the hospital. 

Anthropometry 

Participants' weight was measured using the electronic 
scale produced by Weight-Tronix (Scale Electronics De- 
velopment, New York, USA) to the nearest 0.1 kg. Height 
of the participants was measured in their bare feet using a 
stadiometer (Holtain, Crosswell, Wales, UK) to the nearest 
0.5 cm. Body mass index (BMI) was defined as weight in 
kilograms divided by height in meters squared. Waist cir- 
cumference (WC) was measured at the level of the belly 
button and hip circumference (HC) was measured at the 
largest circumference point around the hips and butt. 
Both were measured using a standard tape measure to the 
nearest 0.1 cm. Each anthropometric measurement was 
performed by a single well trained observer in the same 
room. All subjects, wearing a hospital cotton/polyester 
blend robe with minimal underwear, were measured in 
private. Waist circumference was measured to the nearest 
centimeter at the level of the umbilicus with both arms 
hanging freely at the end of gende expiration, hip circum- 
ference was measured at the spina iliaca anterior superior. 



By the same observers, the circumferences measurement 
variability was within 1 mm with tapes calibrated weekly. 

Skinfold measurement 

Skinfold thickness was measured using the Lange skinfold 
caliper (Cambridge Scientifc, Cambridge, MD) at ten dif- 
ferent locations on the body, including the abdomen (ab- 
dominal skinfold, AS), the supraUiac (suprailiac skinfold, 
SS) [16], armpit (armpit skinfold, ARS), the triceps (triceps 
skinfold, TS), the back (back skinfold, BS), the subscapula 
(subscapula skinfold, SUS), the thigh (thigh skinfold, 
THS), the chest (chest skinfold, CS), the leg (calf skinfold, 
CAS) and the chin (chin skinfold, CHS). Each site was 
measured three times and the average value for each site 
was used. These measures were taken on the left side of 
the body early in the morning to standardize the condi- 
tions of fluid balance [17]. 

Computed tomography 

We used a 64-slice CT Scanner (Somatorn Sensation 64 
CT system), (Siemens Corp., Germany) together with the 
operating software (Software Version syngo CT2005A) to 
scan the abdominal area. Each participant laid down in 
the center of the CT scanning platform and the lumbar re- 
gion was scanned. Scanning was performed with voltage 
of 120Kv, a tube current of 120mAs and an x-ray beam 
width of 1.5 mm; scanning time was 0.5 s. Additionally, 
slice thickness was 5 mm, and images were obtained at 
2 mm intervals. The image reconstruction kernel index 
was B20. The image was furthered processed using the 
commercial software 3D-Doctor Ver. 3.5 (Able software 
Corp., USA). Procedures and specifications recommended 
by previous literature were used to scan images at the 
navel height and to color the areas of abdominal visceral 
fat and abdominal subcutaneous fat followed by an esti- 
mate of the fat areas of the coloring areas (i.e. visceral fat 
area, hereafter denoted as VFAct) [5] . The threshold range 
for adipose tissue in CT was approximately (-260 ± 3) to 
(-10 ± 3) Hu. Before the analysis, CT scans were performed 
to L4-L5 abdominal visceral fat area two times with a three 
day interval on five participants to determine the reliability 
of the measurements. 

Experiment procedure 

The experiment began 2:00 pm every afternoon during 
June 2008 to May 2011. The measuring sequences were 
body weight, height, waist circumference, hip circumfer- 
ence, skinfold thickness, and CT scan. The CT, anthropo- 
metric, and skinfold measurement were all performed by 
the same radiologist and research assistant. 

Statistical analysis 

The demographic statistics of visceral fat area (VFA) 
were measured by CT (VFA^x) are presented using 
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mean ± SD for age, anthropometric measures (i.e., 
height, weight, BMI, waist and hip circumferences) and 
skinfold measurement in the modeling and validation 
groups by CT measurements. The minimum and max- 
imum values were reported in parentheses. 

Multivariate stepwise linear regression was used for 
age, sex, and other anthropometric predictor variables 
(height, weight, BMI, WC, HC, WHR, AS, SS, ARS, TS, 
BS, SS, THS, CS, CAS, CHS) and VFAct was used as 
the response variable. Stepwise regression analysis with 
forward (Fin = 4.00) and backward (Pout =3.99) independ- 
ent variable selection was used for all of the anthropo- 
metric variables in the modeling group (MG) to obtain 
the estimated variables. The regression coefficient, 
standard estimate error (SEE), coefficient of determin- 
ation was further used to construct the VFA estimate 
equation (VFAant)- When independent variables are 
closely related to one another, we may consider remov- 
ing the variable with variance inflation factor (VIF) > 3 
from the estimate equation. Bland-Altman Plot [18] was 
performed to evaluate the goodness-of-fit between the 
results of VFAant estimate equation and VFAct in the 
modeling and validation group (VG). Further, pre- and 
post-test CT data from five patients were used to con- 
firm the test-retest reliability of the measurements. Stat- 
istical significance was set at P < 0.05 for all tests. All 
analyses were carried out using SPSS for WINDOWS 
(Version 16.0; SPSS Inc, Chicago). 

Results 

The subjects were randomly divided into the MG (2/3 total 
subjects) and VG (1/3 total subjects). The modeling group 
consist 87 males (age: 28.56 ±11.20 years, BMI: 24.68 ± 
4.04 kg/m^) and 65 females (age: 33.28 ± 16.23 years, BMI: 
23.15 ± 3.93 kg/ui); the validation group consisted of 43 
males (age: 26.13 ± 12.63 years, BMI: 25.27 ± 3.02 kg/m^) 
and 32 females (age: 27.07 ±8.11 years, BMI: 24.00 ± 
5.60 kg/ra). The results of the VFAct in male (51.49 ± 
40.71 cm^) and female (58.80 ± 37.88 cm^) participants 
were normally distributed. The VFAct and other circum- 
ference measurement results are indicated in Table 1. The 
correlation coefficient and Cronbach's alpha of the CT 
scans measured from the five participants' VFAct were 
0.99 and 0.99 respectively. 

The Correlation coefficient between VFAct and pre- 
dictor variables were as follows: Age (r = 0.70), Height 
(r=-0.13). Weight (r = 0.44), BMI (r = 0.63), WC (r = 
0.75), HC (r = 0.54), WHR(r = 0.74), AS (r = 0.80), SS 
(r = 0.72), ARS (r = 0.81), TS (r = 0.57), SUS (r = 0.76), BS 
(r = 0.76), THS (r = 0.49), CAS (r = 0.45), CHS (r = 0.66), 
CS (r = 0.74), Sex (r = -0.06). ARS and Sex are the highest 
and lowest predictor variables in VFAct (Table 2). 

Age, height, weight, BMI, Sex, WC, HC, WHR, and 
subcutaneous fat thickness (all sites) were selected as 



predictor variables and VFA as dependent variable for 
multivariate stepwise linear regression analysis. Age, 
WC, and AS were the three predictor variables that 
composed the best model for predicting the VFA. The 
estimate equation is shown in Eq. (1). 

VFAant = -144.66 + 1.84 Xi + 1.35 + 0.52 X3 

(r = 0.92, SEE = 14.58 cm^, P < 0.001, n = 152) 

(1) 

The variance inflation factor (VIF) of WC, AS, and 
Age as predictors for VFAant were 1.6, 2.6, 1.9 respect- 
ively. All VIF values for WC, Age, and AS were less than 
3 indicating that there was no bias due to collinearity. 
Assuming a medium effect size, 3 predictor variables, 90 
observations (100 preferred) were needed to reach 80% 
power at a 0.05 significance level [19]. After applying Eq. 
(1) to residual analysis (not shown in text), no significant 
trends were reported. Other than the three predictor 
variables, we added another three variables (BMI, WHR, 
BS) acquired through regression analysis by gradually 
adding the variables into linear regression analysis. After 
performing the analysis six times, a better understanding 
of the collinearity and coefficient of determination be- 
tween variables were observed in Table 3. When BMI, 
WHR or BS were added as additional predictor variables 
in the present model, there was a slight increase in the 
correlation coefficient for VFAant and a slight decrease 
in the SEE. VIF increased from 2.6 to 8.4 which demon- 
strated high collinearity. As result, the predictor variables 
that demonstrated collinearity are not to be included into 
the estimate equation. 

Figure la shows a scatter plot of VFAant estimates 
and VFAct values of the modeling group. The regression 
line is VFAant = 0.83 VFAct + 7.95, r = 0.92, P < 0.001. 
Figure lb shows the agreement in VFA between VFAant 
and VFAct by using the Bland and Altman analysis in 
modeling group. The mean differences in VFA between 
the VFAant and VFAct methods were ± 28.80 cm^. 

We applied Eq. (1) to compute the VFAant for the 
validation group, and compared the estimates with 
VFAct values (Figure 2b). The regression line was 
VFAant = 0.80 VFAct + 10.89, r = 0.92, P < 0.001. Figure 2b 
represented the agreement in VFA between VFAant 
and VFAct and was compared by using the Bland and 
Altman Plot analysis in the validation group. The mean 
differences in VFA between the VFAant and VFAct 
methods were -0.37 ± 33.06 cm^, respectively. 

Discussion 

In the present study, the modeling group was sufficient 
(« = 152) for the proposed analysis. VFAant estimator 
(r > 0.92, P < 0.001) also fared well when compared with 
prior studies (e.g., Naqaiet et al. [20] (r = 0.90-0.92), 
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Table 1 Physical characteristics of the study participants^ 







Modeling group (n = 152) 






Variable 


Male 


(n = 87) 


Female 


(n = 65) 




Age (year) 


28.56± 11.20 


(/S.7, 65.5) 


33.28 ± 16.23 


(18.5, 74.8) 


« 


Height (cm) 


173.39 ±7.85 


(165.0, 195.5) 


160.84 ±6.24 


(163.6,174.0) 




Weight (kg) 


74.41 ±14.17 


(50.0, 121.6) 


59.66 ±9.35 


(52.0, 98.0) 


*** 


BMI (l<g/m^) 


24.68 ± 4.04 


(16.3, 39.9) 


23.1 5 ±3.93 


(16.2, 37.8) 


* 


Waist circumference (cm) 


8241 ±10.82 


(69.0, 122.5) 


81.11 ±10.39 


(68.0, 115.0) 


NS 


Hip circumference (cm) 


98.50 ±7.98 


(87.5, 120.0) 


97.81 ±7.86 


(87.5, 125.5) 


NS 


ACSAcT (cm) 


478.61 ±142.62 


(284.4,1004.8) 


44848 ± 1 22.94 


(288.7,875.2) 


NS 


VFAct (cm) 


5 1.49 ±40.71 


(6.5, 1 79.5) 


58.80 ±37.88 


(13.0, 218.8) 


NS 


SFAcT (cm) 


11 2.56 ±97.66 


(2.7, 513.2) 


165.59 ±99.25 


(20.5, 488.8) 




Waist-hip ratio (WHR) 


0.83 ± 0.06 


(0.72, 1 .02) 


0.83 ± 0.06 


(0.73, 1 .08) 


NS 


Subcutaneous fat thicl<ness (SFT) 












Abdominal, AS (mm) 


18.34 ±13.29 


(2.5, 59.5) 


23.91 ±14.12 


(1.0, 62.0) 




Suprailiac, SS (mm) 


18.56 ±13.33 


(1.8, 59.8) 


24.02 ± 14.11 


(1.0, 62.5) 




Armpit, ARS (mm) 


1 5.67 ± 1 1 .48 


(4.0, 59.0) 


23.88 ± 14.37 


(3.0, 55.0) 


.« 


Triceps, TS (mm) 


16.08 ±11. 77 


(4.0, 58.0) 


24.07 ± 14.51 


(3.0, 56.0) 


*** 


Subscapular, SUS (mm) 


15.99± 11.53 


(4.0, 60.0) 


23.72 ± 14.58 


(3.0, 58.0) 


*** 


Back, BS (mm) 


15.91 ±11.56 


(4.0, 59.0) 


24.02 ± 14.53 


(3.0, 56.3) 


*** 


Thigh, THS (mm) 


1241 ±8.63 


(3.0, 38.0) 


1 9.63 ± 14.30 


(3.0, 77.0) 


*X* 


Calf, CAS (mm) 


1 2.65 ± 8.84 


(3.0, 42.0) 


1 9.07 ± 1 2.92 


(3.5, 55.0) 


»,» 


Chin, CHS (mm) 


12.69 ±8.77 


(2.5, 40.5) 


1 9.02 ± 1 2.97 


(3.0, 56.0) 


»,» 


Chest, CS (mm) 


12.58 ±8.72 


(2.8, 40.2) 


1 9.25 ± 12.72 


(3.2, 55.0) 


». 






Validation group (n = 75) 








Male 


(n = 43) 


Female 


(n = 32) 




Age (year) 


26.13 ±12.63 


(18.2, 70.8) 


27.07 ±8.11 


(18.9, 46.4) 


NS 


Height (cm) 


174.64 ±8.51 


(159.4, 197.0) 


161.35 ±5.53 


(156.0,173.5) 


*** 


Weight (kg) 


77.03 ±10.1 9 


(61.5, 102.0) 


62.58 ± 15.52 


(47.6, 1 06.0) 




BMI (kg/m^) 


25.27 ±3.02 


(20.3, 32.6) 


24.00 ± 5.60 


(17.3, 38.0) 


NS 


Waist circumference (cm) 


84.82 ±10.00 


(68.0, 1 09.0) 


82.36 ± 13.00 


(67.6, 1 22.0) 


NS 


Hip circumference (cm) 


100.1 3 ±5.68 


(88.0, 113.6) 


99.64 ± 1 1 .64 


(86.0, 129.0) 


NS 


ACSAct (cm) 


485.40 ± 1 1 9.58 


(323.7, 795.6) 


47240 ± 1 64.72 


(294.7,995.2) 


NS 


VFAct (cm) 


53.09 ± 40.84 


(8.2, 1 66.2) 


58.98 ±36.1 7 


(18.7, 171.1) 


NS 


SFAcT (cm) 


123.26 ±90.85 


(10.5, 388.0) 


171.34± 118.66 


(8.1, 506.8) 


NS 


Waist-hip ratio (WHR) 


0.85 ± 0.06 


(0.73, 1 .00) 


0.82 ± 0.06 


(0.73, 0.96) 


NS 


Subcutaneous fat thickness (SFT) 












Abdominal, AS (mm) 


18.33 ±14.64 


(3.0, 57.0) 


26.02 ± 14.82 


(3.0, 60.0) 




Suprailiac, SS (mm) 


18.53 ±14.56 


(3.0, 56.5) 


26.42 ± 14.76 


(3.0, 55.0) 




Armpit, ARS (mm) 


17.90 ±14.27 


(3.5, 61.0) 


24.58 ± 13.92 


(3.0, 55.0) 




Triceps, TS (mm) 


18.27 ±14.24 


(4.0, 60.5) 


24.39 ± 13.72 


(2.0, 55.0) 


NS 


Subscapular, SUS (mm) 


18.24 ±14.12 


(3.0, 62.0) 


24.47 ± 13.81 


(4.0, 57.0) 


NS 


Back BS (mm) 


18.14±14.19 


(3.5, 61.2) 


2448 ± 1 3.77 


(3.0, 55.7) 




Thigh, THS (mm) 


12.80 ±9.37 


(2.5, 42.0) 


20.11 ±11.23 


(2.0, 52.5) 
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Table 1 Physical characteristics of the study participants^ (Continued) 

Calf, CAS (mm) 1 3.08 ± 9.54 (3.0, 42.2) 2041 ± 1 2.54 (2.5, 52.0) 

Chin, CHS (mm) 12.89 ±9.1 7 (3.0,41.0) 20.09 ± 12.43 (3.0,53.0) 

Chest, CS (mm) 12.92 ±9.34 (2.8,41.8) 20.20 ±12.38 (2.5,52.5) 

^All values are x ± SD; minimum and maximum in parentheses; 

^Significantly different from male (one-factor ANOVA): NS, not statistically significant; 

*P < 0.05; **P < 0.0 1 ; ***P < 0.001 ; 

BMI; Body mass index; SFA, subcutaneous fat area; VFA, visceral fat area; ACSA, abdominal cross-sectional area. 



Demura et al. [21] (r = 0.87- 0.90), Ryo et al. [22] (r = 
0.88)). Our estimator also had a lower SEE (14.58 cm^) 
when compared with Brundavani et al. [15] (38.7 and 
29.0 cm^ for males and females respectively). The re- 
sults of Bland & Altman plot from the VFA estimate 
equation in MG did not show any systematic error. In 
the VG, the standard deviation (-0.37 cm ) was near 
zero with data sets not normally distributed, which indi- 
cates neither systematic nor proportional error of the 
model was observed. 

Several factors affect how a regression model performs 
[23]: (a) cross-validation, (b) a model is deemed unreli- 
able if its VIF is greater than 10, (c) the sample size, (d) 
the choice of independent predictor variables, (e) and 
the prediction error or standard error of estimate (SEE). 
Anthropometric measures were used because they were 
easy to collect and also biologically relevant as inde- 
pendent predictor variables for estimating VFAant- 



Previous studies focus on VFA estimation used an- 
thropometric measures such as WC, HC, WHR, sub- 
cutaneous fat (AS, CS), age, and BMI separately or 
collectively. While VFA estimation was easy and fast by 
using a single predictor variable, this lead to a higher es- 
timation error. On the other hand, VFA estimation is in- 
convenient when the model used too many predictor 
variables, or relied on expensive technologies. For ex- 
ample, the VFA estimator developed by Demura and 
Sato [21] used total torso fat that can only be obtained 
using Dual-energy X-ray absorptiometry (DXA) or Bio- 
electrical impedance analysis (BIA). Furthermore, total 
torso fat measurement using BIA was difficult [24]; its 
large measurement error in turn also affects the accur- 
acy of VFA estimation. Estimating torso fat weight using 
BIA was a difficult task [24] so there was even larger 
error when used in VFA estimation. Demura et al. [21] 
also proposed using sagittal diameter, WC, subcutaneous 



Table 2 Correlation coefficient matrix between VFAct and predictor variables in the modeling group (n = 152] 





VFAcT 


ARS' 


AS' 


sus' 


BS' 


WC 


CS' 


WHR 


SS' 


Age 


CHS' 


BMI 


TS' 


HC THS' 


CAS' 


Wt Sex 


ARS' 


.81" 
































AS^ 


.80" 


.86" 






























SUS' 


.76" 


.82" 


.78" 




























BS^ 


.76" 


.83" 


.80" 


.92" 


























WC 


.75" 


.63" 


.56" 


.57" 


.58" 
























cs^ 


.74" 


.84" 


.80" 


.80" 


.74" 


.44" 






















WHR 


.74" 


.64" 


.58" 


.59" 


.55" 


.80" 


.54" 




















SS' 


.72" 


.76" 


.83" 


.73" 


.74" 


.48" 


.74" 


.48" 


















Age 


.70" 


.60" 


.67" 


.62" 


.56" 


.28" 


.79" 


.42" 


.68" 
















CHS' 


.66" 


.76" 


.74" 


.73" 


.73" 


.48" 


.72" 


.50" 


.72" 


.54" 














BMI 


.63" 


.57" 


46" 


.50" 


.54" 


.87" 


.34" 


.55" 


.37" 


.16" 


43" 












TS' 


.57" 


.68" 


.69" 


.66" 


.59" 


.31" 


.74" 


.37" 


.56" 


.37" 


.16 


.43*" 










HC 


.54" 


.44" 


.38" 


.38" 


.44" 


.87" 


.22" 


.39" 


.34" 


.09 


.31" 


.86" 


.16" 








THS' 


.49" 


.59" 


.65" 


.61" 


.56" 


.25" 


.69" 


.34" 


.69" 


.64" 


.66" 


.11 


.77" 


.10 






CAS' 


.45" 


.54" 


.61" 


.55" 


.50" 


.20" 


.67" 


.28" 


.65" 


.62" 


.61" 


.10 


.74" 


.07 .85" 






Wt 


.44" 


.38" 


.22" 


.35" 


.35" 


.77" 


.15 


.44" 


.12 


-.02 


.18* 


.84" 


.09 


.81" -.10 


-.08 




Sex 


-.06 


-.06 


-.23" 


.01 


-.07 


.06 


-.07 


.10 


-31" 


-.14 


-.26" 


.15 


-.11 


.01 -37" 


-.22" 


.48" 


Ht 


-.13 


-.18" 


-.30" 


-.13 


-.18* 


.12 


-.24" 


-.05 


-.35" 


-.27" 


-.34" 


.06 


-.23" 


.23" -.38" 


-.33" 


.58" .66" 



^subcutaneous fat thickness; AS: abdominal skinfold, SS: suprailiac skinfold, ARS: armpit skinfold, TS: triceps skinfold, BS: back skinfold, SUS: subscapula skinfold, 
THS: thigh skinfold, CS: chest skinfold, CAS: calf skinfold, CHS: chin skinfold; WC: waist circumference; HC: hip circumference; WHR: waist-to-hip ration; 
Sex(0: female, 1: male); Ht, height; Wt, weight; BMI, body mass index; *, **.Correlation is significant at the 0.05, 0.01 level (2-tailed). 
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Table 3 Multiple regression analysis results for waist circumference (WC) measured with anthropometric measures as 
predictor variable and VFAct as response variable (Modeling group) 



Cumulative dependent variables used in model (n= 152) 



Waist 


Age 




AS^ 


BMI 


WHR BS^ 


Intercept 




SEE (cm^) 


r^ 


2.72 ±0.21 (1.0)" 












-1 67.84 ± 


1 7.33" 


27.296 


0.528 


2.14±0.13(1.ir 


1.64 ± 


0.10(1.1)" 








-1 65.98 ± 


10.36" 


16.322 


0.833 


1.84 ±0.1 5(1 .6)" 


1.35 ± 


0.13(1.9)" 


0.52 ±0.1 5(2.6)" 






-144.66 ± 


1 1 .72" 


15.741 


0.844 


1.61 ±0.27(5.1)" 


1.38 ± 


0.13(1.9)" 


0.52 ±0.1 5(2.6).. 


0.67 ± 0.67(4.4) 




-1 42.54 ± 


11.91" 


15.741 


0.846 


1.1 2 ±0.35(8.4)* 


1.33 ± 


0.13(2.0)" 


0.47 ±0.1 5(2.6)" 


1.10 ±0.68(4.8) 


88.75 ±37.89(3.4)* - 


-1 82.27 ± 


20.96" 


15.517 


0.852 


1.1 2 ±0.35(8.4)' 


1.33 ± 


0.13(2.0)" 


041 ±0.23(6.5) 


1 .08 ± 0.69(4.8) 


88.30 ± 38.03(3.4)* 0.08 ± 0.23(5.5) 


-181.51 ± 


21.14" 


15.565 


0.852 



Regression coefficient estimate ± SEE (VIE, variance inflation factor); r^, determination coefficient; *P < 0.05; **P < 0.001; ' subcutaneous fat thicl<ness; AS: abdominal 
skinfold, BS: back skinfold. 



fat thickness (CS, AS, SUS) and left leg and torso fat 
fractions as predictor variables for VFA estimation. Simi- 
lar to the present model, their approach also used WS 
and AS. Nonetheless, left leg and torso fat fractions were 
difficult to measure [24]. Their model also had a higher 
number of predictor variables. These issues seriously 
limit the practical applications of the VFA estimation 
model. 



Adding these additional predictor variables made the 
model more complex and there was only modest gain in 
the accuracy of VFAant estimation. Thus, we recom- 
mend using only WC, Age, and AS for VFAant estima- 
tion. Bonora et al. [14] also pointed out that when more 
than three predictor variables were used for VFA estima- 
tion, there was no obvious decrease in SEE. For some 
scenarios, it is possible to increase SEE. The adjusted r 



160 



120 



o 80 



< 

^ 40 




VFAj^T= 0.825 VFAct + 7.948 
n=152 
r= 0.92 
P<0.001 



40 



80 120 
VFAcT (cm^) 



160 



200 



40 

E 

3, 20 

H 

O 

< 

I- 

z 

< -20 
> 

-40 



o« o 



4 




o 



o 

• •• 

o 



Male 
Female 



O 



40 



80 120 160 
VFAcT (cm^) 

Figure 1 Modeling group (a) Scatter chart and regression 
line (b) Bland-Altman Plot (mean: 0.0 cm^ mean - 2 
SD: -28.80 cm^ mean + 2 SD: 28.80 cm^). 



200 



200 
160 



E 120 
o 



< 

> 



80 
40 
0 




VFA,„T= 0.804 VFAcT + 10.893 
n=75 
r=0.92 
P<0.001 



50 



100 150 
VFAcT(cm2) 



200 



250 




100 150 
VFAcT{cm2) 

Figure 2 Validation group (a) Scatter chart and regression 
line (b) Bland-Altman Plot (mean: -0.37 cm^ mean - 2 
SD: -33.44 cm^ mean + 2 SD: 32.70 cm^). 
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decreased gradually or approached a fixed constant. We 
observed the same effect in the present study. 

Among the anthropometric predictor variables exam- 
ined in the study, many predictor variables (ARS, AS, 
SUS, BS, WC, CS, WHR, SS, Age) have high correlation 
(r > 0.7) with VFAqx- Other studies have confirmed an 
existing high positive correlation between VFAct and 
WC and WHR. These predictor variables have been 
widely used in VFA estimation studies. Due to the high 
correlation between WC and WHR, the two predictor 
variables are rarely used in the same estimation model, 
lest coUinearity affects the regression estimate. 

Age was used for VFA estimation in Seidell et al. [10], 
Naqai et al. [20], Demura and Sato [21]. These studies 
determined correlations between VFA and age for male 
and female participants were 0.54 and 0.62 respectively, 
which are similar to the data presented in the study. 
These three predictor variables can also be obtained eas- 
ily. In summary, VFA estimation based on Age, WC, and 
AS is highly accurate, reproducible, and practical. 

Several published studies on VFA estimate equations 
also consider sex as one of the predictor variables [15,21]. 
However, in the present study, the correlation between 
sex and VFAct showed a weak correlation between partic- 
ipants (r = -0.06). Based on the results of the correlation 
coefficient matrix, weight was the only highly correlated 
independent variable when compared with sex (r = 0.48). 
The study further confirmed the effect of sex in VFA esti- 
mate by dividing male and female participants 'independ- 
ent variables and acquiring male and female VFA estimate 
equation (not shown in text) via regression analysis. When 
comparing male and female data sets (together and apart), 
the predictor variables, determination coefficient, regres- 
sion coefficient, and SEE values were fairly close and no 
significant differences were found. Hence, sex was omitted 
from VFA estimate. 

There was a high positive correlation (r > 0.7, P<0.01) 
between VFAct and subcutaneous fat thickness at the AS, 
ARS, SUS, BS, CS and SS sites. These six predictor vari- 
ables were also highly correlated with each other (r > 0.7, 
P < 0.01). Thus, multivariate stepwise regression analysis 
was chosen to retain AS and exclude the other five pre- 
dictor variables. 

The study has several limitations. First, the partici- 
pants in the study were young (mean age: 29.2 ± 
12.9 years) and the proportion of VFA >100 cm^ was 
small (14.5%), the present study may have difficulty 
adapting the regression equation for middle-aged people 
who have a higher proportion of VFA >100 cm^ than 
young people. Second, the data are limited to the Tai- 
wanese population which may have different VFA char- 
acteristics than other populations. Third, the number of 
participants was another limiting factor to the study. 
During the recruitment process, we interviewed each 



participant by making sure they are all well and healthy 
before recruiting them into the study, thus only allowing 
these participants may be the cause of the decreased 
sample size. Fourth, the subjects in the study were well 
and healthy individuals and cannot be applied to the 
clinical environment. The main focus of the study is to 
apply the simplified VFA measuring method to healthy 
individuals and contribute to the field of preventive 
medicine. In the future, extending the research to more 
participants and different age groups may enhance the 
understanding of the research. 

Conclusions 

This study demonstrated that rapid and accurate VFA 
estimation can be achieved using only Age, WC, and AS. 
Our approach provides an easy and reliable estimate that 
can be applied widely in health and epidemiology. Using 
too many anthropometric predictor variables will not in- 
crease the precision for VFA estimation. 
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